Spectral properties of acoustic black hole radiation: Broadening the horizon 
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The sensitivity of the black hole spectrum when introducing short distance dispersion is studied in 
the context of atomic Bose condensates. By considering flows characterized by several length scales, 
we show that, while the spectrum remains remarkably Planckian, the temperature is no longer fixed 
by the surface gravity. Rather it is determined by the average of the flow gradient across the horizon 
over an interval fixed by the healing length and the surface gravity, as if the horizon were broadened. 
This remains valid as long as the flow does not induce nonadiabatic effects that produce oscillations 
or some parametric amplification of the flux. 
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I. INTRODUCTION 

In the hydrodynamic approximation, i.e. for long 
wavelengths, the propagation of sound waves in a moving 
fluid that crosses the speed of sound is analogous to that 
of light in a black hole metric [1]. However, Hawking 
radiation relies on short wavelength modes [2-4] . There- 
fore the dispersion of sound waves must be taken into 
account when computing the spectrum emitted by an 
acoustic black hole. To this end, Unruh wrote a disper- 
sive wave equation in a supersonic flow [5]. Through a 
numerical analysis, he then showed that the spectrum 
was robust, provided the dispersive scale ^, the "healing 
length" , is much smaller than the surface gravity scale 
1/k which fixes the Hawking temperature Th = n/2Tr (in 
units where h ~ ^ c ^ 1). This insensitivity was 
then confirmed by analytical [6-9] and numerical [10-13] 
methods. However it turned out that the deviations with 
respect to the standard flux are much harder to charac- 
terize, and so far there is no consensus on what are the 
relevant parameters. 

To address this question, we consider velocity profiles 
characterized by several scales, and numerically compute 
the spectrum by varying them separately. Whenever 
Wmax/Te ^ 10, where cjmax is a critical frequency [12, 13] 
that scales with 1/^ but also depends on the asymptotic 
flow velocity, the flux remains Planckian to a high accu- 
racy, even when the temperature completely differs from 
Th- In fact, we show that the temperature is determined 
by the average of the velocity gradient across the hori- 
zon over a critical length fixed by f and k. In [14], we 
analytically explain the origin of that length. 

These properties hold for large classes of flows. They 
cease to be valid only if the flow induces nonadiabatic 



effects that interfere with the Hawking effect, thereby in- 
troducing oscillations in the spectrum. This is illustrated 
by considering undulations similar to those appearing in 
white holes' flows [15, 16]. For definiteness we use the 
Bogoliubov-de Gennes (BdG) equation in the context of 
atomic Bose condensates. However our results apply to 
other analogue models, e.g. with subluminal dispersion 
relations, and, more generally, to all dispersive theories. 

We have organized the paper as follows. We first 
present the multiparameter fiow profiles we shall use. 
We then analyze the phonon spectrum obtained by nu- 
merically solving the mode equation. In Appendixes we 
briefiy present the concepts which are needed to obtain 
the phonon flux, starting from the BdG equation. 



II. CHOICE OF FLOW PROFILES 

In this paper we consider elongated condensates that 
are stationary flowing along the longitudinal direction x. 
We also assume that the transverse dimensions are small 
enough that relevant phonon excitations are longitudinal. 
In that case, one effectively deals with a 1+1 dimensional 
field. Then, using the analogy [1, 17] between sound 
propagation in the hydrodynamic regime and light, the 
flow defines a metric of the form 



ds'' = -c^dt'^ + [dx - V dtf 



(1) 
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where v is the fiow velocity and c the speed of sound. For 
stationary fiows, c and v only depend on x. Assuming 
that the condensate flows from right to left (w < 0), a 
sonic horizon is present where w = c + v crosses 0. It 
corresponds to a black hole horizon when d^w > 0, and 
a white one when dxW < 0. We shall call it a Killing 
horizon since the norm of the Killing field dt vanishes 
when w = [17]. Its location is taken at a; = 0. 

When ignoring short distance dispersion, the spectrum 
of the upstream phonons spontaneously emitted from a 
black hole horizon is simple, and strictly corresponds to 



2 



the Hawking radiation [1 . 5] . It follows a Planck law and 
the temperature is fixed by the decay rate ^ 



(2) 



For relativistic fields indeed, kk gives the late time de- 
cay rate of the frequency of outgoing modes, as seen 
by asymptotic observers, see Eq. (32.8) in [19] for clas- 
sical waves, and Eq. (3.44) in [4] for Hawking radia- 
tion. As a result, kk fixes the temperature to be exactly 
7h = 'tK/27r, in units where 7i = fc^ = 1, but indepen- 
dently of Co, the value of the sound speed at the horizon. 

When including dispersion, these results are no longer 
true because the total amount of redshift saturates. Nev- 
ertheless, the KK-decay law is recovered near the horizon 
where w = c+v is linear in x. Moreover, when the healing 
length is much smaller than cq /kk , this guarantees that 
the spectrum and the correlations of the Hawking pairs 
are not significantly affected by dispersion [6, 20, 21]. 

As already mentioned, the difficulties start when one 
tries to characterize the deviations from the Planck spec- 
trum at temperature Th. This is because these have dif- 
ferent origin and are governed by different parameters. 
In spite of this, we shall see that they can be analyzed 
and understood. To this end, we shall work with flow 
profiles governed by several scales, and of the form 



c{x) + v{x) = w{x) = wb{x) + Wp{x), 



(3) 



where wb is a reference background profile and Wp a per- 
turbation of smaller amplitude: < |z/;p(a;)| ^ |u)B(a;)|. 
The metric (1) is then completely fixed by introducing q, 



c{x) = Co + (1 - q)w{x), 
v{x) = ~Co + qw{x), 



(4) 



which specifies how c-t-i; is shared between c and v. In all 
simulations, we work with q = 1/2, because it minimizes 
the scattering between left and right moving phonons 
which is not related to the Hawking effect (see Fig. 11 
in [13]). Hence working with q ~ 1/2 will ease the iden- 
tification of the other, more intrinsic, deviations with re- 
spect to the standard flux. Similarly, to have well-defined 
asymptotic modes, we work with infinite condensates and 
with c, V possessing well-defined values for x — > ±oo. 
The background profile is given by 



wb{x) 



Co 



D X sign(a;) 



tanh 



/ kb\x\ 
V coD 



(5) 



The quantity D fixes the asymptotic value of wb , whereas 
n governs the smoothness of the transition from the re- 
gion where wb is linear in x to the asymptotic regime 



^ To adopt the standard terminology, we shall call kk the sur- 
face gravity [18], even though the latter has a dimension of an 
acceleration, and has no clear physical meaning for sound waves. 
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FIG. 1: Upper panel: background flow wb (solid line) for 
n — 3, D = 0.5 and two perturbations Wp (amplified by 
a factor of 5, dashed and dotted lines), with Up = 1, and 
{Hip/ KB, Dp) equal to (0.4, 0.02) and (-0.4, 0.01). Central and 
lower panel: first and second derivative of wb (solid line) and 
of u; = Wb + Wp (dashed and dotted lines) for the same Wp. 
Narrower perturbations induce higher modifications of deriva- 
tives of w. 



where it is constant. We shall work with n = 1 or 2 to 
have a well-defined range ~ cqD/kb where wb is linear, 
and yet avoid the nonadiabatic effects [10, 12] when the 
transition is too sharp, i.e., when n > 4. 

We shall consider two types of perturbations, sym- 
metric ones with respect to the horizon, i.e., Wp{x) 



-x), and then, in Sec. IV, more general ones with- 
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out that symmetry. The symmetric ones are similar to 
WB^■ 



Wp{x) 



Co 



Dp X sign(a;) 



tanh 



Co Dp 



l/rip 



(6) 



The condition < \wp\ < \wb\ is simply < Dp < Z). In 
Fig. 1, we plot w{x), its first and second derivatives for 
different values of k^/kb and Dp. 

The smoothness parameter Up will range from 0.5 to 
4. We shall see that it only affects the flux marginally. 
In fact the most important quantity is Dp. Its role is to 
fix 



Fiii, 



Co Dp 



(7) 



the interval of x over which w{x) is linear. In this respect, 
it should be noticed that the surface gravity is equal to 

KK = 9x(c + W)U=0 = KB + Kp, (8) 

irrespectively of the value of Dp, and that of rip. For 
relativistic fields, Dp would play no role since the tem- 
perature is KK/27r (when ignoring gray body factors, a 
correct approximation in our settings when q = 1/2). 
For dispersive fields instead. Dp plays a crucial role. 
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FIG. 2: log]^g(/ij/rH) versus uj/kb for different values of Dp, 
but the same surface gravity kk . The solid line represents the 
Planck flux at temperature Tn = kk/^tt. When Dp = 0.1, 
the phonon flux hardly differs from it. For smaller values 



of Dp, the flux is weaker. 



and for Dp — >■ 0, it is flxed by 



the background surface gravity kb- The flxed parameters are 
A/kb = 15, D = 0.4, Kp/hiB = 0.5, n = Up = I. 
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III. SPECTRAL ANALYSIS 

In a stationary condensate flowing with a velocity v, 
the dispersion relation of phonons is given by [22] 

[CO - vkf = n\k) = c^e + ^ = c^e + 4? , (9) 

4m^ 

where w is the conserved frequency dt, rn the atom mass, 
and A the frequency associated with the healing length 
^ = h/V2mc by A = V^c/C 

In what follows we study the properties of , the mean 
occupation number of phonons of frequency uj sponta- 
neously emitted to the right region, i.e. against the flow, 
when the initial state is vacuum. In order to focus on new 
results, we have decided to present in the Appendixes a 
summary of the concepts needed to compute n^^ , starting 
from the BdG equation applied to the flows of Eq. (3). 
The figures of this paper have been obtained by numeri- 
cally solving Eq. (B6) with the code of Ref. [13]. 

Our program can now be clearly defined: we shall 
study in the perturbed flows for increasing Dp to 
see how the phonon field "responds" to the introduc- 
tion of Wp. We start with the symmetrical profiles of 
Eq. (6), and then consider asymmetrical ones. The spec- 
trum in the unperturbed flow of Eq. (5) has been studied 
in [10, 12, 13]. 



In Fig. 2 wc present the energy flux 

fuj^ujn^, (10) 

as a function of w for different values of the amplitude 
Dp, keeping fixed all the other parameters. When work- 
ing with a relativistic massless field, one would obtain a 
Planck spectrum with a temperature Th — 'tK/27r, irre- 
spectively of the value of Dp. To ease the comparison 
with the dispersive fiuxes, we have represented this fiux 
by a solid line. Working with the BdG equation (B6), 
we see instead that the fiux varies with Dp. In fact it 
increases monotonically. Moreover, when Dp — > 0, we 
verified that it coincides with the fiux one would obtain 
in the reference flow wb- We also see that when Dp is 
sufficiently large, the flux saturates and agrees to a high 
accuracy with the Planck spectrum at temperature Th, 
in accordance with the robustness [5, 6, 10] of the flux. 

Since comparing energy fluxes f^j is not convenient, we 
shall use the temperature function defined by 

1 , , 

naj = , I . -. 11) 

exp(w/T<^) - 1 

As such, is simply another way to express How- 
ever it presents the great advantage of being constant 
whenever the flux is Planckian. 

In fact, as can be seen in Fig. 3, for all values of Dp, 
remains remarkably constant, until lo approaches the 
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FIG. 3: Tui/Tb of Eq. (11) versus uj/kb for various Dp but 
the same value of kk = 1.5kb. Tb ~ kb/Stt is the Hawking 
temperature of the background flow. The values of the pa- 
rameters are those of Fig. 2. For all values of Dp, T^j hardly 
varies until u — >■ Wmax of Eq. (A3), where the flux vanishes. 
Since ojmax ^ 10 Tb, as can be seen from the horizontal upper 
scale, these spectra are, to a high precision, Planckian. 

critical frequency Wmax of Eq. (A3). At that frequency, 
the signal drops down as it must do, since the Bogoli- 
ubov transformation (B8) becomes trivial, because the 
negative norm mode ip-u no longer exists above Wmax- 
This cutoff effect relies on the asymptotic value of the 
flow V + and has been studied in detail in [12, 13]. In 
what follows we focus instead on the deviations of the 
spectrum which depend on the near horizon properties 
of u + c. 

To quantify the small deviations from a Planck spec- 
trum at temperature 

To = lim T^, (12) 

we studied the relative deviation 

A^.;) ^ (13) 
Jo 

evaluated for u = Tq. For all values of Dp we found 

|AT(ro)| < 3. X 10-4. (14) 

We note that At (To) 7^ is not due to numerical errors, 
it really characterizes the deviations. 

Equation (14) implies that, to a high accuracy, the 
spectrum remains Planckian even when Tq completely 
differs from Th- (In the present case, it differs by 33%, 
since kk = 1.5 kb.) The smallness of Ay (To) also im- 
plies that the parameters that govern the deviations from 
Planckianity differ from those governing the shift Tq— Th- 



FIG. 4: To/Tb versus Dp for Kp/kb from 0.2 to 1, and with 
A/kb — 15, D = 0.3, n — Up — 1. Along each curve the sur- 
face gravity kk of Eq. (8) is constant. It fixes the temperature 
only for D 3> D"'*. The vertical line gives the value of the 
healing length ^. 



Hence when studying the robustness of the Hawking 
flux against introducing dispersion, one must differen- 
tiate these two types of deviations. 

These are our flrst important results. They generalize 
what was observed in [12, 13] to a much wider class of 
flows, namely, first, that the thermality is well preserved 
whenever Wniax/Ts ^ 10, and, second, that the relative 
temperature shift Tq/Th — 1 is much larger than At{Tq). 
In the sequel we work with Wmax/Te ?J 10 satisfied, and 
study how the low frequency temperature Tq depends on 
Kp, Dp, Hp and A. 



B. The critical value of Dp 

In Fig. 3, we saw that Tq starts from K'q/2tt for Dp — s- 0, 
and asymptotes to Th ~ K}^/2t: for Dp larger than some 
critical value (We have checked that this is also 

the case when Kp < 0, and this down to Kp = — kb, in 
which case the surface gravity kk vanishes.) To under- 
stand what fixes i'p"', we first consider series of fiows 
wb + Wp{Dp) for different values of Kp at fixed kb- For 
each series, the surface gravity kk = kb + Kp is thus con- 
stant. The resulting temperatures Tq{Dp) are presented 
in Fig. 4. This figure reveals many interesting features. 
For small values of Dp, Tq increases linearly with Dp in 
a manner essentially independent of Kp. In addition, for 
each series, when Dp is large enough, Tq saturates at the 
standard result Th = kk/Stt. 

Given the linear behavior of To for small Dp and the 
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FIG. 5: To/Tb versus Dp for A/kb from 12 to 20, and with 
Kp/kb = 0.5, D — 0.3, n — Hp = 1. The vertical lines give the 
values of the healing length ^ = \/2co/A. One can see that 
Dp"* of Eq. (17) is not (directly) related to 



saturation for large values, i.e., 

2ttTo{Dp < Df) ^ KB 
27rTo(I?p » Df) ^ KB 



sDr, 



(15) 



wc define Z?"'* as the value of Dp where these intercept, 
and obtain 



(16) 



To understand what fixes the slope s, we consider series 
of flows with different values of the ultraviolet disper- 
sion scale A. The results are presented in Fig. 5. It is 
numerically easy to see that s oc h?/^ . Hence we get 



(17) 



We now present a simple argument telling how D^^'^ 
should scale with the various parameters. In [13], by 
considering wb of Eq. (5), it was observed that the lead- 
ing deviations due to dispersion are governed by inverse 
powers of Wmax/'^B- For these unperturbed flows, Wmax of 
Eq. (A3), even though initially defined as the frequency 
such that the turning point xt-p{ijj) of Eq. (A2) is rejected 
at — oo, also corresponds to the frequency such that the 
turning point is at the edge of the domain ^ cqD/kb, 
where wb is linear in x. 

When working with w = wb + Wp of Eq. (6) and 
with Dp <^ D, the linear domain is now given by 
of Eq. (7). Hence we expect that the deviations will be 
now governed by inverse powers of wmp/'^k, where WMp 
is the new critical frequency fixed by xiin. The value of 



WMp is such that the turning point |a;tp(w)| of Eq. (A2) 
equals xiin. Using Eq. (7), la;'P(a;M,^)l = xnn gives 



i^Mp = A ( —D. 

fip 



''P(^Mp)| 
3/2 



Xlin gives 



(18) 



Since the deviations from the standard result should be 
small when ujmp/kk ^ 1, I?"'' should scale as 



d: 



crit 



CX 



(kK A2)1/3- 



(19) 



A rigorous study [14] of the connection formulas [6-9] 
confirms the result of this simple reasoning. 

Equation (19) is in agreement with Eq. (17) for Kp <C 
Kb . In the following section, we shall numerically validate 
it for all values of Kp. Before proceeding, we checked that 
Dp '* does not significantly depend on the background 
quantities n, D and kb, as one could have expected. 
The weak, subleading, dependence on Up is studied in 
Sec. HID. 



C. The averaged surface gravity 

So far we have some understanding of the Planckian 
character of the fiux in the robust regime when Dp > 
Dp"*-, and of the value of Dp"* [14]. Instead, what fixes 
outside this regime, for Dp < Dp"*, is terra incognita. 

We shall now numerically establish that, to a good 
approximation when Wmax/TB ?i 10, the temperature is 
determined by the average of the gradient dw/dx over a 
width that we parametrize by d^ = d^ + d^ , where d^ 
and d^ are respectively the width on the left and on the 
right calculated from the Killing horizon at a; = 0. 

In fact, given that the dispersion of Eq. (9) defines a 
'healing length' ^ = a/2c/A, it is not surprising that ^ 
in turns defines a minimal resolution length, such that 
details of w{x) smaller than d^ are not "seen" by the 
phonon field. In other words, because of dispersion, it is 
as if we were dealing with a horizon of width d^.^ More- 
over, when d^ ^ d^ , this means that the center of the 
effective horizon is displaced with respect to the Killing 
horizon to Xhor = {df — d^)/2. 

To test the idea that the temperature be determined 
by an average "surface gravity," we introduce 



— / dx 



dw{x) _ w{df)-w{~dl) 



H J-d 



dx 



(20) 



^ This important result should be opposed to the possibility dis- 
cussed in [9, 23] according to which T^j could be governed by 
a local, (xi-dependent, function of w, e.g. by the value of dxW 
evaluated at the turning point of Eq. (A2). In fact our numeri- 
cal results indicate that it is a non-local quantity that fixes the 
temperature, and in a manner essentially independent of u). This 
is in sympathy with the ideas that quantum gravitational effects 
might blur the horizon [24, 25] and that dispersion might regulate 
the black hole entanglement entropy [26] . 
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FIG. 6: Dots: Tq/Tb versus Dp computed for A/kb = 15, 
D = 0.4, Up — n — 1, Kp/fvB = 0.5. Solind line: best fit curve 
obtained by using R,/2Tr of Eq. (21). The perfect agreement es- 
tablishes that the temperature To equals R,/2n. Dashed lines: 



asymptotic behaviors for Dp Dp'^' 
Eqs. (15) and (16)]. 



Using Eqs. (5) and (6), we get 

Wp{df)+Wp{df) 
t = Kb H ; 



and Dp > D^"' [i 



O 



Den 



2n 



where we assume that <C D/kb, so that the average 
of the background term wb is approximately kb- 

We proceed as follows: We compute To for about 200 
values of Dp for some fixed values of the perturbed flow 
parameters {up, Kp) and the dispersive scale A, i.e., along 
series as those represented in Figs. 4 and 5. Out of 
this set we extract five quantities, namely, kb^^'\ np'^*^*'\ 



Kp^^^\ df, and d?, by fitting k(kb, 



,d^,d^;Dp)/2n 



of Eq. (21) with a nonlinear least-squares method. The 
agreement between the couples Tq{Dp) and the fitted 
function R{Dp)/2iT is really striking — see Fig. 6 for an 
example. 

Moreover, the fitted values kb*-*^*-* and Kp^*^*) are in 
very good agreement (less than 1.5% and 1%, respec- 
tively) with their values used in the flow w. This is a 
necessary check to validate that it is the average of dxW 
that governs the temperature. However, the agreement 
between rip^**^ and Up is less good. This indicates that 
the actual temperature is not exactly given by the aver- 
age of Eq. (20). ^ Even though this is a subdominant 



effect, it is interesting to see under which conditions a 
better agreement between Up^^'^'^ and Up is reached. 



D. Subleading effects in Up 



The parameter Up of Eq. (6) determines the smooth- 
ness of the transition between the region where Wp{x) is 
linear, with constant derivative Kp, and the region where 
it is fiat, equal to DpCo- Larger values of Up make that 
transition sharper. Therefore higher derivatives of w will 
become very large at the end of the linear-regime region, 
near xnn of Eq. (7). As a consequence, nonadiabatic ef- 
fects become more important, thereby generating larger 
deviations between Tq and R/2tt of Eq. (21). 

To investigate these deviations, we compare Tq{Dp) for 
various values of Hp. As it appears from Fig. 7, left panel, 
s and do not significantly depend on Up. However, 
the transition region between the linear regime and the 
constant regime depends on Up. The transition is indeed 
monotonic for small values of rip whereas, for rip > 1, 
To oscillates around Th = kk/^tt. The monotonic func- 
tion of Eq. (20) will of course never describe these os- 

(fit) 



cillations. However, for rip < 1, one can compare rip 
with rip. The result is shown in Fig. 7, right panel. The 
dots and the squares are computed using different values 
of A/kb; respectively 15 and 20. One sees that rip^' is 
always larger than rip. One also sees that the discrep- 
ancy becomes smaller when rip decreases, as expected, 
since nonadiabatic effects are smaller for smoother pro- 
(21) files, hence smaller rip. One should not overestimate the 
relevance of this discrepancy because we found that the 
value of rip*^*^*) is correlated in the fitting procedure to the 
difference d? — dp. For instance, when imposing d? 



(fit) 



d^ 



is reduced by a factor of about 2. 



The important conclusion is the following. When 
nonadiabatic effects are small, the temperature of the 
phonon fiux is well approximated by l/27r times the av- 
erage of dxW over a certain width d^ [see Eq. (20)]. 



E. The scaling laws of the width 



Unlike kb , 



tp. Dp which characterize the flow w and 



A which fixes the healing length, d^ is not a parameter 
entering the mode equation (B6). Therefore the value 
obtained in a fit cannot be compared with some a pri- 
ori known value 
emergent, effective, length scale 

Before using numerical results to determine how it 
scales when changing kk, Hp and A, we point out that 



In fact d^ should be conceived as an 



^ We notice that different averaging procedures could have been 
considered in the place of Eq. (20), and they would have given 



slightly different results. At this point, we are not in the position 
to distinguish them. What is important is that the scaling laws 
we shall derive in Sec. Ill E are independent of the particular 
choice. 
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FIG. 7: Left panel: To/Tb versus Dp for different values of Up of Eq. (6) with Kp/kb = 0.5, D — 0.4, n = 2. To monotonically 
approaches Th for Up < 1. When Up > 1 instead, it exceeds Th near Dp"^'' and then approaches it with oscillations. Right 



panel: the relative difference {up 



(fit) 



with respect to A/kb (dots, A/kb = 15; squares, A/kb = 20). 



Up)/ Up versus Up. The difference decreases for smaller rip, in a rather insensitive way 



we can predict how it should do, as is closely related 
to D"/'*- of Eq. (19). Indeed, taking the Hmits of Eq. (21) 
for small and large Dp we obtain 



h{Up <^ ) ^ Kb ^ 3 — , 

Co 

R{Dp > - KK, 

Co 



and hence, using Eq. (15) and (16), we get 



(22) 



df = —. 



(23) 



In other words, the critical value of Dp separating the two 
regimes of Tq{Dp) shown in Fig. 6 strictly corresponds to 
the width one should use in Eq. (20) to follow To(Z)p). 
Using Eq. (19), we get 



df — d X 



Co 



(^kA2)i/3' 



(24) 



by a linear ^ least-square fit of 



log(d{— ) = log(rf) + alog( — 
Co V '*B 



&log( ^ 



clog 1 



Kb 



(26) 



We used a grid of 110 values of (k^. A). The correspond- 
ing values of d^^Hp, A) are extracted from series of about 
200 couples {Tq,Dp) using the procedure described after 
Eq. (20). Kp/KB ranges from 0.1 to 1, with step 0.1 and 
A/kb ranges from 10 to 20 with unitary step. The other 
parameters are D = 0.4, n — Up — 1, q = 0.5.^ 
The result of the fit is 



a = -0.627, 
c = -0.33, 



b = 0.024, 
d = 1.58, 



(27) 



where the numerical error is on the last figure. The ex- 
ponent a differs from about 5% from its expected value 



where d is a constant. In terms of the healing length 
C = y2co/A, this gives 



-1/3 



(25) 



where dx = cq/kk is the surface gravity length. Hence, 
when ^ dx, is significantly larger than ^. 

The prediction of Eq. (24) can be numerically validated 



* When the function is linear in the parameters, a close- form so- 
lution of the minimization problem is available. This method 
is therefore preferable to non-linear ones. 

^ Because of numerical noise, for some series (A,Kp), it was not 
possible to obtain a sufficiently large number of points (To, Dp). 
In particular it is difficult to obtain Tq for large values of Dp. 
We therefore decided to use only the series for which we obtain 
To for every value of Dp from to 0.1 (namely 100 points, with 
step 0.001). In this way, there are points with Dp > D"'' in the 
series (To, Dp) used to determine d^. This condition is satisfied 
by 71 series over 110. 
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—0.667. The exponent b is also close to its expected value 
0. Finally, c perfectly matches its theoretical estimate. 
Furthermore, we notice that the unknown constant d in- 
troduced in Eq. (24) is of the order of 1. 

In brief, for the considered class of flows, we have es- 
tablished that, to some accuracy, 

• the spectrum is Planckian, 

• the temperature equals K/27r for a width dj, and 

• scales according to Eq. (24) with d ~ 1.6. 

These are the main results of this paper. 

From our fits we also found that the broadened horizon 
is slightly displaced with respect to the Killing horizon 
by a shift |a;hor| ~ O.ldj. However, because of the sym- 
metry of Eq. (21), we cannot know the sign of ihor- To 
determine it, it is necessary to consider asymmetric flows. 



IV. ASYMMETRIC FLOW PROFILES 
A. Monotonic perturbations 

In this section, we consider profiles of the form 

Wps{x) = Wp{x - 5) - Wp{~5), (28) 

which are simply Wp of Eq. (6) shifted from the horizon 
by 6. We have subtracted the constant term so that 
wb + Wps still vanishes at a; = 0. In Fig. 8, Tq is plotted 
as a function of Dp for different 5 (left plot), and as a 
function of i5 for different Dp (right plot). The behavior 
of these curves is quite complicated. To understand it, 
we first consider smooth perturbations, i.e. perturbations 
with Dp > D;"' « 0.06. 

From the right panel of Fig. 8, when Dp > Dp"^, but 
for both signs of S, we see that To monotonically interpo- 
lates from Tfj = KK/27r, when S is small and the perturba- 
tion is close to the horizon, to Tb = kb/'Ztt, when 6 — > oo. 
In that case, the perturbation is far from the horizon and 
thus no longer contributes to kk which equals its back- 
ground value Kb ■ When considering instead well localized 
perturbations, that is values of Dp smaller than D"'*, we 
see that Tq oscillates, but only for 6 sufficiently negative, 
i.e. when the perturbation is sufficiently displaced in the 
region where the decaying mode oscillates. ^ The origin 
of these oscillations can be understood as follows. The 
scattering on the well localized bump in dxw{x) interferes 



with the scattering on the horizon, thereby producing os- 
cillatory behaviors. 

It is quite clear that the average of Eq. (20) cannot de- 
scribe these oscillations. However this should not be con- 
ceived as limiting the validity of the observation that the 
temperature is given by this average. Indeed, these os- 
cillations result from interferences between two spatially 
separated scatterings. In these circumstances it is mean- 
ingless to apply Eq. (20), which is based on the near 
horizon properties of w. 

In addition, when oscillations are found for Tq as a 
function of Dp or (5, the Planck character of the spectrum 
is also lost. Indeed, the parameter At(ix') of Eq. (13) 
characterizing deviations from Planckianity also displays 
oscillations (of relative amplitude ^ 20%) (see Fig. 9). 
As can be seen, the number of peaks increases with in- 
creasing values of —5. This is to be expected since more 
"room" is available to fit resonating modes. In fact this 
phenomenon can be viewed as a precursor of the black 
hole laser effect [28]. If the perturbation were strong 
enough to behave as a white horizon, the resonating 
modes would start growing exponentially [29, 30]. In that 
case complex-frequency modes appear when the distance 
between the horizons is sufficiently large, and their num- 
ber increases with that distance. 

To complete this analysis of the flows of Eq. (28), we 
verified that R of Eq. (20) still approximately governs 
the temperature. In Fig. 10, both the numerical data 
(dots) and the fitted function (solid line) are reported 
for the series Tq{5), with Dp = 0.05 and 0.15, of Fig. 8 
(right panel) . This function is obtained as a two- variable 
fimction of 5 and Dp by using the numerical data of To- 
The fitting region is restricted to those values of 5 for 
which there are no oscillations and for three values of 
Dp, namely Dp = 0.05,0.1,0.15. As a reference, the dot- 
ted line represents KK((5)/27r, i.e. the temperature one 
would obtain without dispersion, for Dp = 0.05. As a 
check, we also fitted the values of kb, Hp and Up and 
compared them with the corresponding constants used 
in the numerical analysis. The relative difference is very 
small, respectively, around 0.1%, 1% and 5%. 

As one can see in Fig. 10, the agreement between the 
actual temperature Tq and the average surface gravity 
R/2'k is quite good in the parameter region where there 
are no interference effects. One can also sec that for 
Dp = 0.05, R/2-K of Eq. (20) follows Tq much more closely 
than the surface gravity KK('5)/27r. 

In summary, the analysis of asymmetrical perturba- 
tions first indicates that the shift ihor is always negative. 
That is, the broadened horizon is displaced towards the 



We have adopted this formulation because, for subluminal dis- 
persion, we expect that the sign of x^^r wiU be the opposite, 
and this in virtue of the 'symmetry' between the behavior of 
modes for sub and super-luminal dispersion [14]. Indeed, for 
(sub) super-luminal dispersion the turning point is in the (sub) 
super-sonic region, see Eq. (A2). In both cases, the mode decays 
on the other side of the turning point [10, 12], 



This is very similar to what has been found in inflation. When in- 
troducing a sharp modification of mode propagation before hori- 
zon exit, there are interferences between this localized scattering 
and the standard mode amplification which produce oscillations 
in the power spectrum, see [27] for references and a critical anal- 
ysis on the generic character of this. 
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FIG. 8: On the left panel, Tq/Tb as a function of Dp for various values of S, and, on the right, Tq/Tb as a function of S with 
Dp from Dp = 0.02 to Dp — 0.3. A/kb = 15, D = 0.4, ftp/ftB = 0.5, n = Up = 1. On the left panel, for values of <5 > —0.05, 
To{Dp) behaves as in Fig. 6, whereas, for more negative values of S, To behaves differently. On the right panel, the asymmetry 
of To under S ^ -6 is clearly visible for Dp < Dp"^ ^ 0.06. 



region where the decaying mode oscillates. Moreover we 
find that the value of the shift does not depend directly 
on 6. By this we mean that Xhoi- depends on 5 only be- 
cause of Eq. (24) is a function of 6 through kk{S). 
Second, for sufficiently smooth perturbations we found 
that the temperature is still approximately given by the 
average of Eq. (20). Finally, when the perturbation is 
sharp and sufficiently deep inside the supersonic region, 
we found that the spectrum is affected by interferences, 
and explained why it is so. 



B. Undulations 

Supersonic flows possessing a sonic horizon with a 
negative surface gravity kk behave like the time re- 
verse of black hole horizons: long wavelength modes are 
blueshifted when scattered on such a horizon. Thanks to 
dispersion, the blueshift effect saturates, and as a result, 
the scattering matrix is well-defined. Moreover, for each 
cj, this matrix is closely related to Eq. (B8) evaluated 
in the corresponding black hole flow [12, 13]. In spite of 
this correspondence, it has been recently observed that 
white hole flows display a specific phenomenon related to 
the limit w — >■ [15, 16]. Indeed, these flows produce a 
zero-frequency mode, a stationary undulation, which has 
a finite wavelength kz and a high amplitude. 

To schematically describe such undulation and study 
its impact on the emitted spectrum, we consider the fol- 
lowing perturbations: 

^^^^Axsm{kzx-9)xe~^^/^^\ (29) 

Co 



A is the amplitude of the undulation, and 9 fixes its phase 
on the horizon at a; = 0. The new length A is taken to 
be larger than l/kz, and comparable to Dcq/kb which 
characterizes the linear regime of the background flow 
wb- 

In Fig. 11; To is plotted versus 9 for different values 
of kzi keeping fixed Akz so that the surface gravity 
of Eq. (8) does not change with kz- Indeed one gets 
hk{9) = kb + Akz cos{9). When k^ ^ ^/d^i as expected 
from our previous results, the temperature is very close to 
kb/^tt, that of the background flow wb- Hence, indepen- 
dently of the value of 9, the perturbation is completely 
washed out by the dispersion. On the contrary, when the 
wavelength is larger than the width of the horizon, i.e. 
when kz < 27r/c?^ w SOkb/cq, we expect that 2ttTo will 
be approximately given by kk(^) given above. This can 
be observed for kz equal to Skb/cq. One notices that 
the curve is not exactly symmetric around 9 = tt. But 
this is also expected from the former analysis where we 
saw that moving the perturbation in the supersonic or in 
the subsonic region affects differently the spectrum. 

The novelties are found for intermediate values of kz 
(say, between IOkb/cq and 20kb/co). For these values, 
6z ~ 2TT/kz is comparable to d^. Hence the undula- 
tion can enter into resonance with the scattering on the 
horizon, thereby producing a kind of parametric ampli- 
fication of the flux. This can be observed for 9 = 
and kz = 10kb/co where the temperature is larger than 
KK/27r = 1.5Tb. The fact that this resonance occurs for 
wavelengths ~ c?^, reinforces the fact that the effective 
width of the horizon is indeed d^. 
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FIG. 9: The parameter At(i^) of Eq. (13) versus Lojn-s, for 
different values of &. For 5 = 0, continuous line, the perturba- 
tion is symmetric with respect to the horizon, and At(i^) is 
extremely small in agreement with Eq. (14). When 5 > 
and (—5) < rf{ ~ 0.2co/fi:B, AT(a;) stays very small and 
no oscillations are found. Instead when —J ^ dj, oscilla- 
tions of large amplitudes appear even for low frequencies, i.e. 
oj <^ tJmax — 13 Tb, thereby indicating that the spectrum is 
no longer Planckian. The fixed parameters are A/kb = 15, 
D = 0.4, D-p = 0.02, Kp/^B = 0.5, n = = 1. 



FIG. 11: To/Tb as a function of Q for Eq. (29), varying kz 
and A so that Akz = 0.5kb stays fixed. A = 2Dco/k,b, 
A/kb = 15, D = 0.4, 71 = 1. For very high values of kz, i.e. 
kz = 100 and 30 kb/cq, the temperature is approximately 
k,b/'2'ji, i.e. equal to R/2ti of Eq. (20). For small values of kz, 
i.e. kz = 5kb/co, To also approximately follows R(6)/2ti ~ 
HK{d)/2Ti since the undulation wavelength is larger than d^. 
In the intermediate regime instead, for kz = 10 and 20kb/co, 
the temperature widely oscillates and no longer follows k/2tt. 
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FIG. 10: Numerical data (dots) for the series To{S)/Tb with 
Dp = 0.15 « 2.57?™* and 0.05 « 0.8D"" of Fig. 8, right plot. 
The solid lines are the values of the fitted function of Eq. (20). 
The dotted line represents the surface gravity k,k{5)/2h for 
Dp = 0.05. The agreement between the temperature and 
R/2ti is quite good in the domains with no interference. The 
lowering of the peak and its shift to the left for small values 
of Dp are rather well reproduced. None of these effects are 
found when using the surface gravity kk(5). 



V. CONCLUSIONS 

By considering series of supersonic flows consisting on 
local perturbations Wp defined on the top of a background 
flow wb , we were able to identify the parameters that flx 
the spectral properties of Bogoliubov phonons emitted 
by a (black hole) sonic horizon. 

We show that the flux remains Planckian to a high ac- 
curacy, see Eq. (14), even when the temperature strongly 
differs from Th ~ the standard value fixed by the 

surface gravity, sec Fig. 3. This result implies two things. 
In spite of the fact that each mode of frequency oj has 
its own turning point of Eq. (A2), the Planckianity first 
implies that the temperature function of Eq. (11) is 
not determined by some function of wix) evaluated at 
a w-dependent location, as Tj^ is common to all low fre- 
quencies with respect to Wmax of Eq. (A3). Second, it 
implies that the parameters governing the small devia- 
tions from Planckianity differ from those governing the 
difference between the low frequency temperature Tq and 
Th. In fact, whereas the difference Tq — Th is governed 
by the near horizon properties, deviations from Planck- 
ianity are governed by inverse powers of Wmax/To, where 
f^max is fixed by T^as = T* + Dp, cliaracterizing the global 
properties of the supersonic transition region, such as the 
asymptotic velocity excess. 

Moreover, in contradistinction to the relativistic case 
where the temperature is locally fixed by kk of Eq. (8), 



11 



we show that the temperature Tq is determined, to a high 
accuracy when Wmax/Tn ^ 10, by a nonlocal quantity: 
the spatial average across the horizon of dxW over a width 
dj, as if the horizon were broadened. ® Hence, whenever 
w{x) hardly varies over in the near horizon region, S of 
Eq. (20) equals the surface gravity kk, and the standard 
temperature is found. Instead, when w varies on scales 
shorter than d^, Tq is no longer the standard one, but it is 
well approximated by K/27r for all frequencies uj <C i^max- 
This result is reinforced by the fact that the effective 
horizon width scales according to a well-defined law 
given in Eq. (24). 

Furthermore, the above results are basically valid for 
all flows. They cease to be valid under the same restricted 
conditions, namely when the flow contains a sharp per- 
turbation well localized far enough from the horizon. 
Then the perturbation induces a scattering that inter- 
feres with the Hawking effect, thereby engendering os- 
cillations. These are similar to those found in inflation- 
ary spectra, and are also attributable to interferences be- 
tween two scatterings. It is clear that the local average 
of Eq. (20) cannot describe these oscillations which arc 
nonlocal on a scale much larger than the width d^. 

Hence we can conclude that when Wmax/TH ^ 10 and 
when the scattering essentially arises near the horizon, 
the spectrum is near Planckian, and with a temperature 
approximately given by k of Eq. (20). 

Finally, to mimic what is found in white hole flows [15, 
16], we studied profiles containing an undulation on the 
top of a smooth profile. By varying its wavelength, 
we have found three different behaviors that match the 
above analysis. 
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Appendix A: Asymptotic modes and turning points 

We briefly show why and how the number of modes of 
frequency uj depends on the asymptotic properties ofw = 
c + V. This is important as it explains why, when there 



® This is in agreement with what was found in [29] when con- 
sidering the black hole laser effect [28]. In that case, the effect 
disappears when the distance Su between the black and the white 
horizon becomes too small, because the dispersive field is unable 
to 'resolve' the two horizons. When using linear profiles for both 
horizons, one finds that the critical value of (5h — rfj- 



is dispersion, there is no flux above a critical frequency 
'^max whose value is fixed by the asymptotic value w. 

In homogeneous subsonic flows, \v\ < c, for real u>, two 
roots of Eq. (9) are real, and describe the right and 
left moving modes. The other two roots are complex, 
conjugated to each other, and describe modes that are 
asymptotically growing or decaying, say to left. Only 
the first two correspond to perturbations that must be 
quantized, as the last two are not asymptotically bound. 
In infinite condensates, when canonically quantizing a 
field, only asymptotically bound modes (ABIVI) should 
be included in the spectrum [13, 29]. 

In homogeneous supersonic flows, the situation is quite 
different. For real w, there exists a critical frequency 
Wniax below which the four roots are real, as can be seen 
in Fig. 12 (right panel), and above which two arc real and 
two complex, as in subsonic flows. The spectrum of ABM 
thus contains four modes for frequencies < a; < Wmax 
and only two for uj > Wmax- At a; = Wmax the two extra 
roots merge and the curve Q{k) is tangent to the line 
UJ — vk. Hence the group velocity in the lab 

v^, = dku; = dkin + vk), (Al) 

vanishes at that frequency. Solving this equation one 
finds that Wmax is proportional to A, but also depends 
also on the supersonic velocity excess \v\ — c (see Fig. 4 
in [12] for details). 

In inhomogeneous flows that cross once the speed of 
sound and with two asymptotic flat regions, i.e., as those 
of Eq. (5), one gets a new situation: when to > Wmax, 
there are still only two modes, whereas for uj < Wmax, 
there are now three independent ABJVI [12]. There are 
only three independent modes because the semiclassical 
trajectories associated with the two extra roots necessar- 
ily possess a turning point for a; < 0, in the supersonic 
region. Hence these two roots describe the initial and 
final value of the momentum of the same mode. 

The location of the turning point a;*P depends on uj. Its 
value is easily obtains by solving Vg^ = [see Eq. (Al)]. 
Using Eqs. (8) and (9), for cj <C A, one finds [10] 

For — > fixed A, or A — >■ oo fixed uj, the turning point 
coincides with the Killing horizon at a; = 0, as expected 
from the behavior of the semiclassical trajectories in the 
absence of dispersion. [When considering subluminal dis- 
persion by flipping the sign of the fc"* term in Eq. (9), one 
verifies that x*'^{uj) flips sign but keeps the same norm.] 
When UJ — )- cj„iax, the turning point x*p(cj) — >■ — oo. 
Hence, the threshold frequency cjmax is now fixed by the 
asymptotic value of w for x — > — oo. In our velocity fiows 
of Eq. (3), one has 

^max A/(i?as), (A3) 

where Das = -—w{x ~ — oo). Using Eqs. (5) and (6), one 
has Das — D + Dp . 
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FIG. 12: Graphical solution of the dispersion relation (9) for subsonic (left panel) and supersonic flow (right panel). Solid 
curves: ±f2(fc)/K. Left panel, dashed line: {u — vk)/^.. Right panel, dashed, solid, dotted lines: uj — vk, respectively, for 
(jj < CiJniax, 1-0 = cJmax, <^ > i^max ■ For iu < cjmax, two extra (real) roots exist in the left lower quadrant. They cease to exist for 
UJ > cjmax and this explains why there is no Hawking radiation above t^max- 



Appendix B: Mode analysis and scattering matrix 

We present the main concepts that characterize the 
scattering of phonons in atomic Bose condensates in sta- 
tionary flows containing one sonic horizon. We follow the 
treatment of Rcf. [13] where more details can be found. 

Phonon elementary excitations arc described by a 
quantum field operator, (f>. In this work, we shall use 
the relative perturbations defined by 



(Bl) 



where 5' is the atom's field operator, and 'I'o the con- 
densate. The dynamics of (f> is determined by the 
Bogoliubov-de Gennes equation [22]. When using cf> and 
when h ~ 1, one obtains 



idf 



\Tp — ivdx + TOC^] (p + rat? 



(B2) 



where c is the x-dependent speed of sound 

c (x) = , (B3) 



where po{x) = |^'o(x)|^ is the mean density of condensed 
atoms, and g{x) is the effective coupling constant among 
atoms. These functions play no role in the sequel: only 
c{x) and v{x) matter. The quantity 



Tp = --^vdj-dx, 
2m V 



(B4) 



is the kinetic operator that acts on 0. 

In stationary situations, (j) can be expanded in Killing 
frequency eigenmodes 



where the discrete sum over a takes into account the 
number of ABM at each lu. Inserting Eq. (B5) in Eq. (B2) 
yields the system for the couple of modes 4>Z{x)7 'PZi^) 



[-{lj + ivdx) -Tp- mc"] ^Zi^) = mc^cpZix). 



(B6) 



As discussed above, in backgrounds containing one 
sonic horizon, they are two or three ABM depending on 
whether u is larger or smaller than Wmax- Therefore the 
cj component of the field operator reads 



(B7) 



The modes (/)^ and 0^ have positive norm, they are asso- 
ciated with the usual real roots of Eq. (9) (sec Fig. 12), 
and they describe, respectively, propagating right-u and 
Icft-u moving waves with respect to the condensed atoms. 
On the contrary, (ip-^)* has a negative norm. It is asso- 
ciated with the extra roots that exist in supersonic flows 
for Lo < Wmax, and describes phonons that are trapped in- 
side the horizon {x < 0) in the supersonic region. These 
modes describe the negative frequency partners of the 
outgoing Hawking phonons represented by and d^. 

In our infinite condensates, the above modes become 
superpositions of plane waves both in the left and right 
asymptotic regions of Eq. (B7). We can thus construct 
without ambiguity the in- and out-mode bases. The 
in modes are such that each of them contains only one 
asymptotic branch with group velocity directed towards 
the horizon. The definition of out modes is analogous 
and based on the criterion that each out mode contains 
only one asymptotic branch with group velocity directed 
away from the horizon. 

As shown in Appendix A, for w > Wmax, there exist 
only two positive norm modes, and therefore in- and out- 
mode bases are linearly related by a trivial (elastic) 2x2 
transformation. Instead, when tu < Wmax, the in and out 
bases contain three modes which mix with each other in a 
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nontrivial way by a 3 x 3 Bogoliubov transformation [13] 
(^'_"J* = + (^-;^)* + (B8) 

The standard normalization of the modes yields relations 
such as (from the first equation) 

|«^P-|/3-..P + |^|' = l. (B9) 

When the initial state is vacuum, the final occupation 
numbers are given by the Bogoliubov coefficients 

n^ = m\ nl^\B^\\ ^^^^^ 



From these relations one sees that the norm of /3i^ fixes 
the occupation number of outgoing phonons sponta- 
neously produced by the scattering of the quantum field 
near the black hole horizon. In the standard analysis 
without dispersion, n^j is Planckian and describes the 
Hawking effect [4]. In that case, for massless conformally 
invariant fields, n^, the number of left moving phonons 
spontaneously produced, identically vanishes. In the 
presence of dispersion, one generally finds <^ rii^. 
Hence the scattering of a dispersive field on a sonic hori- 
zon basically consists of a more general and slightly mod- 
ified version of the standard Hawking effect. This is 
true when Wmax ^ 2k. For larger values of k, the de- 
viations become large, but the structure of Eq. (B8) and 
the meaning of its coefficients remain unchanged [13, 31]. 
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